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Abstract — In this paper, we propose novel low-complexity 
adaptive channel estimation techniques for mobile wireless chan- 
nels in presence of Rayleigh fading, carrier frequency offsets 
(CFO) and random channel variations. We show that the selective 
partial update of the estimated channel tap-weight vector offers 
a better trade-off between the performance and computational 
complexity, compared to the full update of the estimated channel 
tap-weight vector. We evaluate the mean-square weight error of 
the proposed methods and demonstrate the usefulness of its via 
simulation studies. 

Index Terms — Adaptive filter, channel estimation, carrier fre- 
quency offsets, mean-square weight error, multipath channel, 
Rayleigh fading. 

I. Introduction 

In wireless communications environment, signals suffer 
from multiple reflections while travelling from the transmitter 
to the receiver so that the receiver ends up getting several 
replicas of the transmitted signal. The reflections are received 
with different amplitude and phase distortions, and the overall 
received signal is the combined sum of all the reflections. 
Based on the relative phases of the reflections, the signals 
may add up consttuctively at the receiver. Furthermore, if 
the transmitter is moving with respect to the receiver, these 
desttuctive and constructive interferences will vary with time. 
Therefore, the performance of the wireless systems critically 
dependent on the availability of accurate channel estimation 

m, ®, m. 

Eke, periodic system variations resulting from mismatches 
between the transmitter and receiver carrier generators can be 
damaging to the performance of the estimators, even for small 
carrier frequency offsets 0. Adaptive filtering techniques 
are suitable for tracking of such time variant channels. There 
are many adaptive filter algorithms that are widely used for 
channel estimation, but they either have a high mean-square 
error with slow convergence rate (such as, least mean square 
(LMS) and normalized least mean square (NLMS) algorithms) 
or a high computation complexity with fast convergence rate 
and low mean-square error (such as, recursive least square 
(RLS) algorithm) J6), Q, (8). Thus, mean- square error 
(MSE), convergence rate and computational complexity are 
three important points in selecting of adaptive algorithms for 
channel estimation and this is considered in choice of the 
applied algorithms. To address these problems, Ozeki and 
Omeda [9| published the basic form of an affine projection 
algorithm (APA) and Shin and Sayed developed tracking 



performance analysis of a family of APA in [ 10 1. Furthermore, 
in order to reduce of computational complexity, one possible 
way is that only part of the estimated tap-weight vector is 
selected for update in each time iteration. In ifTTIl . lfl2l and 
lf]~3l adaptive filter algorithms with selective partial update 
method (SPU) has been shown to have low computational 
complexity and good performance. In this paper, based on 
general formalism for the family of affine projection in |10] 
and using the approaches that presented in [11] and lfl2l . we 
present a novel algorithms with general formalism which we 
called it low-complexity family of affine projection algorithms 
based on selective partial update method. Then, we use all of 
these algorithms for estimation of mutipath Rayleigh fading 
non-stationary channels. 

The rest of this paper is organized as follows. Section 
II describe a brief introduction on the statistics of mobile 
wireless channels. In Section III, we will have a breif review 
on some known adaptive estimators and introduce cyclic non- 
stationary channels model. Subsequently, the low-complexity 
family of affine projection algorithms based on selective 
partial update method will be presented in Section IV. Finally, 
Section V analyzes the performance evaluation of the pro- 
posed estimation approaches by computer simulation results to 
demonstrate the effectiveness of the proposed algorithms for 
fading non-stationary channels estimation in mobile wireless 
environments. Section VI concludes this paper. 

II. Mobile Wireless Multipath Channel Model 

The complex baseband representation of the mobile wireless 
channel impulse response can be described by [14J 



h(t, t) = } 7fc(i)x fc (r)(5(T - r k ) 



k=l 



(1) 



where r^. and jj, (t) are respectively the delay and path loss of 
the fc-th path, ccfc(r) is a time-variant complex fading sequence 
of the fc-th path that models the time-variations in the channel. 
Hence, the frequency response at time t is 

r+oo 

H(tJ)± h(t,r)e-^ fT dr (2) 



without loss of generality, the sequence x(t) is assumed 
to have unit variance. Several mathematical models can be 
used to characterize the fading properties of x(r) and con- 
sequently of the channel. A widely used model is known 



as Rayleigh fading. In this case, for each r, the amplitude 
\x(t)\ is assumed to have Rayleigh distribution [3|. Therefore, 
the sequence {xk(r)} are modeled as independent Rayleigh 
fading sequences and the channel is referred to as a multipath 
fading channel. In additional, the atu-correlation function of 
the sequence x(r), now regarded as a random process, is 
modeled as a zeroth-order bessel function of the first kind, 
namely, 

r(k) 4 Ex(t)x(t — k) = J {2nf D T s k). (3) 

This commonly used choice of the auto-correlation function 
is based on the assumption that all scatters are uniformly 
distributed on a circle around the receiver, so that the power 
spectrum of the channel fading gain x(t), in continuous-time, 
would have the following well-known [/-shaped spectrum 

S(J) ^==, l/l < Id (4) 

In equation ([3]), T s is the sampling period of the sequence 
x(t), Jd is called the maximum Doppler frequency of the 
Rayleigh fading channel, and the function Jo is defined by 

Jo = — cos(y sin 9)d9. (5) 

Jo 

The Doppler frequency inversely proportional to the speed of 
light, c = 3 x 10 8 m/s and proportional to the speed of the 
mobile user, v, and to the carrier frequency, f c . 

Due to the motion of the vehicle, 7fc(i)'s are modeled 
to be wide-sense stationary (WSS), narrowband complex 
gaussian process, which are independent for different paths. 
Furthermore, 7fc(i)'s for all k have the same normalized 
time correlation function and different average powers a 2 ,. 
In nonstationary channels, a first-order approximation for the 
variation of a Rayleigh fading coefficient x(t) is to assume 
that x(t) varies according to the auto-regressive model 

x(t) = r(l)x(T - 1) + y/l - |r(l)| 2 77(r) (6) 

where r(l) = Jo(27t/dT s ) and rj(r) denotes a white noise 
process with unit- variance. 

The two-ray 1131 . typical urban (TU), and hilly terrain (HT) 
lfl4l . [ 16 1 models are three commonly used delay profiles. For 
the two-ray profile with equal average power on each ray, the 
delay spread is td/2 (td = T\ — To), i.e., a half of the delay 
different between the two rays. The delay spreads for TU and 
HT delay profiles are 1.06 and 5.04 [is, respectively lUTl . 

III. Overview on Some Known Adaptive Channel 
Estimators 

Considering Fig. 1, x(n) denotes a sequence that is trans- 
mitted over an unknown channel of finite impulse response 
w(n) of order M. d(n) and e(n) are received signal and 
the output estimation error, respectively. Consider the received 
signal d(n) assuming that arisen from the linear model 

d(n) = u(n)w (n) + v(n) (7) 

Wo {n) = w t (n)e^ n (8) 
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Fig. 1. A prototypical adaptive filter setup. 

where wt(n) is an unknown channel tap-weight vector and 
time variant that we wish to estimate it, the term of e J n indi- 
cates the carrier frequency offsets (5), v(n) is the measurement 
noise, assumed to be zero mean, white, gaussian and indepen- 
dent of u(n), and u(n) = [x(n),x(n — 1), x(n — M + 1)] 
denotes 1 x M row input regressor vectors of the channel with 
a positive-definite covariance matrix, B, u = E[u* (n)u(n)]. 

In sequel, we shall adopt a model for the variation of w t (n). 
Many models can be defined for this variation such as random- 
walk or Markovian models [18 |. One particular model that is 
widely used in the adaptive filtering literature is a first order 
random-walk model iflOl . lfT9l . On the other hands, we know 
that the carrier frequency offsets between transmitters and 
receivers and fading phenomenon are two important trails of 
the mobile wireless channels. Therefore, we can use achievable 
model for variations of fading channels until the performance 
of adaptive estimators is enhanced. 

For this reason, the following both cyclic and random non- 
stationary model for variation of channel is introduced Q, 
E) 

w t (n) = w t + E(n) (9) 

where w t is a constant vector that represent the non-fading 
part of the channel, whereas S(n) represent the fading part of 
channel that it can be modeled as autoregressive (AR) process 
of order p [20|. We use common AR(1) model for this process 
as 

E(n + 1) = a£(n) + q(n) (10) 

where 

<| a \< 1 

and following from [8|, the sequence q(n) is assumed to be 
i.i.d, zero-mean with autocorrelation matrix E{q(n)q* (n)} = 
Q and wt (n) have a constant mean which we shall denote by 
w t , Ewt{ri) = w t . In this paper, a and autocorrelation matrix 
of q(n) are Jo(27r/uT s ) and Q = \/T— a 2 1, respectively. 
This non-stationary model of channel is practical in wireless 
communication systems and especially in channel estimation 
and channel equalization applications. 

It is well known that the adaptive update scheme for the 
LMS estimate of w {n) is given by (HI 

w(n + 1) = w(n) + fiu* (n)e(n) (11) 

where 

e(n) = d(n) — u(n)w(n) 



TABLE I 
APA FAMILY ALGORITHMS 



A lu< wit li in 






R 
H 




NLMS 


K = 1 


e = 


p = o 


D = 1 


APA 


K < M 


e = 


13 = 


D = 1 


BNDR-LMS 


K = 2 


e = 


13 = 


D = 1 


R-APA 


K < M 


e^O 


p = o 


D = 1 


PRA 


K < M 


e^O 


/3 = 1 


D = 1 


NLMS-OCF 


K < AI 


e = 


(3 = 


D > 1 



is the output estimation error at time n and is the step size. 

To increase the convergence speed of the LMS estimator, 
the NLMS algorithm was proposed which can be state as [8| 



w(n + 1) = w(n) + 



;U*(n)e(n). 



(12) 



i|u(n)|| 2 

For reducing the computational complicity, the estimated tap- 
weight vector w(n) and the excitation vector u(n) of the 
adaptive estimator can be partitioned into B blocks with length 
L for each blocks that B = Al/L and it shall be an integer 
which are defined as lfl3l 

u(n) = [xi(n),x 2 (n), ...,x B {n)} (13) 

win) = [wf(n),w%(n),...,wl(n)} T '. (14) 

The selective partial update NLMS algorithm for a single block 
update at every iteration can be derived from the solution of 
the following minimization problem IT31 

min \\wj(n + 1) — Wj(n)\\ 2 (15a) 

subject to d(n) — u(n)w(n + 1) (15b) 

by using the method of Lagrange multipliers 12D . the update 
equation for selective partial update NLMS is given by 



where 



Wj(n + 1) = Wj(n) + - — —u*(n)e(n) 



j = argmax||nj(n)|| 



(16) 



for 1 ^ i ^ B and it denotes the number of blocks that should 
be updated at each iteration. 

From [ 1 1 , the general class of affine projection algorithms 
can be stated as 

w(n + 1) = w{n - /3(K - l)) + 

fiU*{n)(eI + U(n)U*(n))- 1 e(n) (17) 

where e(n) — d(n) — U(n)w(n — f3(K — 1)), and 

U(n) = [u(n), u(n - D), u(n - (K - l)D)] T (18a) 
d(n) = [d(n),d(n - D), d(n - (K - l)D)] T . (18b) 

Based on dTvT > and by specific choices of the parameters 
{K, D, e, {3} are resulted in different affine projection algo- 
rithms. The particular choices and their corresponding al- 
gorithms are summarized in Table U For NLMS-OCF, it is 
further assumed that u(n — jD) is orthogonal to u(n), u(n — 



D), ...,u(n — (j — l)D) [ 1 1 . The motivation for using D > 1 
is to increase the separation, and consequently reduce the 
correlation. For PRA, the weight vector is updated once every 
K iterations that K is positive integer and most algorithms 
assume K < M. 

IV. A Low-Complexity Family of Affine Projection 
Algorithms 

This section introduces a low-complexity family of affine 
projection algorithms based on selective partial update method. 
From ll22l . the generic estimated tap-weight vector update 
equation can be stated as 



w(n + 1) = win - P{K - 1)) + nU*(n)C(n)e{n) 



(19) 



where C(n) is the K x K matrix and it is obtained from Table 
UTI This table shows that many classical and modern adaptive 
filter algorithms can be derived through dT9b . 

Considering (12), the constrained optimization problem, 
which is solved by the proposed algorithms are given by 

min \\wh(n + 1) - Wh{n)\\ 2 (20a) 

w h (n+l) 

subject to din) — U(n)w(n + 1) (20b) 

where h = {ti, T2, ts} denotes the blocks that updating at 
every adaptation. By using the lagrange multipliers method 
and considering ( fT9l , the general update equation of low- 
complexity family of affine projection algorithms is given by 

w h {n + 1) = w h (n - K') + fj,U^(n)C h (n)e{n) (21) 

where C\(n) is the K x K matrix and it is obtained from 
Table HI1 K' = j3(K - 1), and we have 

U h {n) = [U Tl {n), U T2 {n), U Ts (n)} (22) 

U TS {n) = [u TS (n),u TS (n-D),...,u TS (n-(K-l)D)} T (23) 



Ui(n) = [ui{n),Ui{n — D),...,iii(n— (K — 1)D)] 



(24) 



In above equations, Uh{n) is the K x LS matrix, U TS (n) is 
K x L, u TS (n) is 1 x L and f7,(n) is the K x L matrix. For 
obtaining the indices of h, we compute the values of 



Ti(Ui(n)U* {n)) for l<i<B 



(25) 



then, S largest blocks are selected for updating. The K x K 
matrix Ui(n)U*(n) is assumed to be full rank (i.e., in- 
vertable). The inequality L > K is a necessary condition 
for Ui(n)U*{n) to be full rank. Note that, for K = 1, 
Ti(Ui(n)U* (n)) reduces to ||ii(n)|| 2 . Thus, d25l l is consistent 
with the selection criterion for the SPU-NLMS algorithm. 

From (I2TI 1. with different choices of the parameters 
{K, D, e, (3} and Ch(n) from Table JTT1 the new adaptive 
estimators such as low-complexity version of BNDR-LMS, 
NLMS-OCF and PRA based on selective partial update 
method will be derived. 



TABLE II 

Correspondence Between Special Cases of (f2TT > and Various Low-Complexity Adaptive Filtering Algorithms 
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C(n) = [Ui^U^n)}- 1 


C h (n) = [U h {n)U* h {n)}^ 


R-APA 


K < M 
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D = 1 


C(n) = [el + U(n)U*(n)]- 1 


C h {n) = [eI + U h {n)Ul{n)]- 1 


PRA 


K < M 


13 = 1 


D = 1 


C(n) = [el + U(n)U*(n)]- 1 


C h {n) = [el + U^U'in)]- 1 
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K < M 


13 = 


D > 1 


C(n) = [Ui^U'in)}- 1 


C h (n) = [^(n)^(n)]- 1 
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Fig. 2. The cumulative distribution function of a Rayleigh fading sequence 
corresponding to a vehicle moving at 80Km/h. 

V. Performance Evaluation by Simulation Studeis 

Extensive computer simulations have been conducted to 
demonstrate the performance of the adaptive estimators for 
multipath Rayleigh fading non-stationary channel estimation. 
In adaptive filtering problems, recursive algorithms are used 
to find the channel tap-weight that produces the minimum 
error between the desired value and received value. Since 
the channel can be modeled as a tap-delay line filter, the 
channel estimation problem can be formulated as that of an 
adaptive filtering problem. Fig. Q] illustrates the process of the 
adaptive filtering and we can use of this process to estimate 
the non-stationary channel taps. The objective of the adaptive 
algorithm is to find the optimal tap-weight vector that gives 
the minimum mean-squared error JS). 

In this paper, we consider a wireless channel with two 
Rayleigh fading rays; both rays are assumed to fade at the 
same Doppler frequency of fo = 10Hz. The channel impulse 
response sequence consists of two zero initial samples (i.e., 
an initial delay of two samples), followed by a Rayleigh 
fading ray, followed by another zero sample, and by a second 



Fig. 3. The first 1000 samples of a learning curve that results from tracking 
a Rayleigh fading non-stationary channel using LMS and family of APA 
algorithms. 

Rayleigh fading ray and then, for another channel response, 
we will add two zero samples after a second Rayleigh fading 
ray. In other words, we are assuming a channel length of M=5 
or 7 taps with only two active Rayleigh fading ray, so that the 
weight vector that we wish to estimate has the form 

[0 x 2 (t) x 4 (t) ] (26) 

or 

[ x 2 {t) x a {t) ] (27) 

Train adaptive algorithm to estimate and track these multipath 
non-stationary channels. Assume random binary phase shift 
keying (BPSK) input signal of unit variance is transmitted 
across the channel and use it to excite the adaptive filter. 
Assume further that the output of the channel is observed in 
the presence of white additive Gaussian noise with variance 
al = 10- 3 , sampling period T s = 0.8^s, a = 0.9, a\ = 1CT 4 
and carrier frequency offsets ^ = 0.0001. Fig. [2] shows the 
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Fig. 7. A plot of a typical tr aject ory of the amplitude of the first Rayleig h ray and its estimate. The plots zooms on the inte rval 1000 < i < 2000. |(a)| 
LMS algorithm with ft = 0.0l[(b)]APA algorithm with K = 1 and fi = 0.l[(c)]APA algorithm with K = 2 and fi = 0.25[(d)]PRA algorithm with K = 3 
and (i = 0.25 1(e)] LC-APA algorithm with K = 3,B = 5,S=1 and fi = 0.25[®]LC-APA algorithm with K = 3,B = 5,S = 2and^t = 0.25 [(g)] 
LC-APA algorithm with if = 3,B = 5,5 , = 3and/x = 0.25. 




Fig. 4. The first 1000 samples of a learning curve that results from 
tracking a Rayleigh fading non-stationary channel using LMS, APA and Low- 
Complexity APA algorithms. 

cumulative distribution function of the amplitude sequence. A 
plot of the learning curve of all adaptive estimators is gener- 
ated by running them for 30000 iterations and averaging over 
100 independent experiments and show the good performance 
of proposed algorithm with lower computational complexity 
that convergence rate of its is comparable with full update APA 
estimators. Figs. [3] and [4] show only the first 1000 iterations 
of a typical learning curves. Fig. [5] shows mean-square error 
(MSE) curves versus step-size for the proposed scheme and 
other algorithms that they are presented in this paper. In this 
simulation, the Doppler frequency is fixed at 10Hz for both 
rays, channel length M = 5 taps, and the MSE values are 
obtained by averaging over 100 independent realizations. 
In Fig. [6] the Doppler frequency vary from 10Hz to 40Hz 



Fig. 5. A plot of the tracking MSE as a function of the step-size for Rayleigh 
fading rays with Doppler frequency at 10Hz. 

in increments of 5Hz and generate a plot of the MSE as a 
function of the Doppler frequency. Run all adaptive estimators 
for 60000 iterations in each case and average the squared- 
error curve over 1000 independent experiments. Also, we used 
the optimum value of the step-size for all algorithms which 
we selected these optimum values from minimizing of mean- 
square weight error function over step-size fi Q. This figure 
shows the robustness of adaptive estimators for different values 
of Doppler frequency. 

Fig. 7 shows typical trajectory of the amplitude of the first 
ray and its estimate by different adaptive estimators. These 

1 In Appendix A, we calculate the mean-square weight error of the proposed 
method over the cyclic non-stationary channel models. For the sake of arriving 
to minimum mean-square error of estimation, we can minimize this function 
over fi. 
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Fig. 6. A plot of the tracking MSE as a function of the Doppler frequency 
of the Rayleigh fading rays 



plots zooms onto the interval [1000, 2000]. The results show 
the good tracking performance of adaptive algorithms for 
multipath Rayleigh fading non-stationary channels tracking 
and proposed method is better than LMS and APA (K = 1) 
algorithms and comparable with APA (K = 3) algorithm but 
our proposed method has lower computational complexity than 
ordinary family of APA algorithms. 

VI. Conclusion 

The tracking performance of the new low-complexity family 
of affine projection algorithms in the presence of carrier 
frequency offsets and random channel nonstationarities was 
studied in this paper. A general formalism for proposed algo- 
rithms was also developed, and the mean-square weight error 
under steady-state conditions was analyzed. The performance 
of our methods were compared with those of some other 
estimators. The proposed schemes presented a good trade-off 
between the performance and the computational complexity. 
The simulation results were included to demonstrate our 
claims. 
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Appendix A 
The Mean-Square Weight Error of the Cyclic 
Non-Stationary Channel Estimation 

To find the mean-square weight error of low-complexity 
family of affine projection algorithms, following from |8|, Q, 
let us praise a definition and a assumption as follows 

Definition: For the cyclic model of channel nonstationarities 
, the weight-error vector is e(n) — Wt(n)e 3 n — Wh(n), and a 
priori and a posteriori estimation errors are defined as e a (n) 
U h {n)e{n-K'), e p {n) = U h {n)(e{n+l)-g(n)e^ n ), respec- 
tively. We also define the error-vector correlation matrix and 



mean-square weight error as Z(n) = E[e(n — K')e* (n — K')], 
T{n) = Tr[Z(n)}, respectively. 

Assumption: The noise v(n) is i.i.d and statistically inde- 
pendent of the regression matrix Uh(n) with variance a 2 = 
E\v(n)\ 2 . Moreover, the random variables d(n), v(n), Uh(n) 
have zero mean and Uh(n)U^(n) is invertible and e(n) is 
statistically independent of Uh(n). Furthermore, The sequence 
q(n) is independent of the initial conditions {wt(— 1), w(— 1)}, 
of the u(n),v(j) for all n,j and of the d(j) for all j < n. 

Thus, by using the above definition and considering © and 
( [Tol l, the weight-error vector recursion can be rewritten as 

£(71 + 1) = e{n-K') + g{n)e^ n -pXJl{n)C h {n)e{n) (28) 

where 

g(n) = (w t + oS(n) + q{n))e^ - (w t + S(n)). (29) 

Introducing the noise vector as 

V T {n) = [v(n) u(n-l) ... v(n - (K - 1)D)]. 

Then, from (O and above assumption, the output estimation 
error can be expressed in term of e a (n) as 



e(n) = e a (n) + V(n). 



(30) 



In continues, considering above assumption and using (l30l ) in 
d28T > and applying it for calculation of error-vector correlation 
matrix Z(n + 1), we get 



Z(n+1) = E[e(n+l)e*(n + l)} 

= Bi(ra) + B 2 (n) + B 3 (n) + B A {n) 



where 

Si(n) 
B 2 {n) 

B 3 (n) 

B 4 {n) 
and 



Z(n) - fj,Z(n)R N + E[e[n - K')g*(n)e- j9n ] 
E[g(n)e^ n e*{n - K 1 )] + E || g{n) || 2 

- liE[g{n)e^ n e*{n- K')] 
-nR N Z(n) + [i 2 Z(n)A(n) 

- f iR N E[e(n - X>*(n) e - J * n ] 
p 2 a 2 v E[U* h (n)C h (n)C* h (n)U h (n)} 



R N = E[U£(n)C h {n)U h (n)} 
A(») =E[Ut{n)C h {n)U h {n)Ut{n)Ct{n)U h {n)]. 

For the sake of obtain the error-vector correlation matrix, we 
require to evaluate the cross correlation between the tap-weight 
error vector and g(n), and E \\ g(n) \\ 2 . At first, by taking 
the squared norm and expectation from both side of d29b , we 
obtain 



E[\\g(n)\\ 2 } =Tr(Q)- 



1 - e 



J* |2 



Tr{W t ) 



Tr{&) 



(31) 



= where 



6 = lim E(E{n)E*{n)). 



From © and O we can find that 

n 



1- a[ 



Then, in steady-state conditions, following the derivation in 
ll23l . we can define that 

E(s(n - K')) = H(n) = He^ n 
E(s(n - K')E*(n)) = Y(n) = YeP n 



where 



H=[I- riR N - e^l^wtil 



Y = [a* {I - (iR N ) ~ e J '*/]- 1 (a*(l - ae^)Q - e jW Q) 



J*/ 



Therefore, using (I29l l and above expressions, we can expand 
B\(n), B2(n) and Bz(n) as following 

Bi(n) = Z(n) - fiZ(n)R N + (e _i * - l)w* t H 

+ {a*e-^ - l)Y 
B 2 (n) = (e J * - l)w t H* + (ae jf - 1)Y* + Tr(Q) 

+ I 1 - e> 9 | 2 Tr(W t )+ I 1 - ae J * | 2 Tr(0) 

- fiR N (e^ - l> t IP - nR N (ae jiS - 1)Y* 
B 3 {n) = -riR N Z(n) + ri 2 Z{n)A{n) 

- fiR N (e- j<!! - l)w*H - LiR N (a*e- j * - 1)Y. 

Thus, by calculating error-vector correlation matrix and apply- 
ing the definition of T(.) in it, we obtain 

T(n + 1) =(1 - 2fxTr(R N ) + fL 2 Tr{A{n)))T{n) 

+ pL 2 a 2 v Tr(E[U* h (n)C h (n)C* h (n)U h (n)}) + F. 

Finally, using steady-state conditions T(/i) = lim T(n, fi) = 

n— too 

lim T(n + mean-square weight error of cyclic non- 

n— >oo 

stationary channel estimation may be given as 

(i 2 * 2 v Tr(E[UZ(n)C h (n)C* h (n)U h (n)}) + F 



T(p) = 
where 



2nTr{R N ) - n 2 Tr(A(n)) 
F =|1 - e^\ 2 Re{Tr(W t (I-2F a ))} 



1 - ae^\ 2 Re{Tr{Q{I-2a*F f} ))} 



Re{Tr(Q(I - 2(a* 



Fp))} 



and 



F a = (/ - fxR N )[I - fiR N - e^I}- 1 
F P = (I- nR N )[a*(I - nR N ) - e^I}- 1 

References 

[1] J. G. Proakis and M. Salehi, Digital Communications, 5th ed. McGraw- 
Hill, 2008. 

[2] T. S. Rappaport, Wireless Communications: Principles and Practice, 
2nd ed. Upper Saddle River, NJ: Prentice-Hall, 2001. 

[3] A. Goldsmith, Wireless Communications. Cambridge Univ Press, 2005. 

[4] M. Rupp, "Lms tracking behavior under periodically changing systems," 
in Proc. of Eur. Signal Process. Conf, Island of Rhodes, Greece, Sep. 
1998. 



[5] A. R. S. Bahai and M. Saraf, "A frequency offsets estimation technique 
for nonstationary channels," in Proc. of ICASSP '97, Munich, Germany, 
Sep. 1997, pp. 21-24. 

[6] B. Widrow and S. D. Stearns, Adaptive Signal Processing. Englewood 
Cliffs, NJ: Prentice-Hall, 1985. 

[7] S. Haykin, Adaptive Filter Theory, 3rd ed. NJ: Prentice-Hall, 1996. 

[8] A. H. Sayed, Adaptive Filters. Wiley, 2008. 

[9] K. Ozeki and T. Umeda, "An adaptive filtering algorithm using an 

orthogonal projection to an affine subspace and its properties," Election. 

Commun. Jpn, vol. 67-A, no. 5, pp. 19-27, 1984. 
[10] H. C. Shin and A. H. sayed, "Mean-square performance of a family 

of affine projection algorithms," IEEE Trans. Signal Process., vol. 52, 

no. 1, pp. 90-102, Jan. 2004. 
[11] T. Aboulnasr and K. Mayyas, "Complexity reduction of the nlms 

algorithm via selective coefficient update," IEEE Trans. Signal Process. , 

vol. 47, no. 5, pp. 1421-1424, May 1999. 
[12] S. Werner, M. L. R. de Campos, and P. S. R. Diniz, "Partial-update nlms 

algorithms with data-selective updating," IEEE Trans. Signal Process., 

vol. 52, no. 4, pp. 938-948, Apr. 2004. 
[13] K. Dogancay and O. Tanrikulu, "Adaptive filtering algorithms with 

selective partial updates," IEEE Trans. Circuits Syst. II., vol. 48, no. 8, 

pp. 762-769, Aug. 2001. 
[14] R. Steele, Mobile Radio Communications. New Yourk: IEEE Press, 

1992. 

[15] Y. Li and N. R. Sollenberger, "Spatial-temporal equalization for is-136 

tdma systems with rapid dispersive fading and co-channel interference," 

IEEE Trans. Veil. Technol, vol. 48, pp. 1182-1194, Jul. 1999. 
[16] S. Ariyavisitakul, N. R. Sollenberger, and L. J. Greenstein, "Tap- se- 
lectable decision feedback equalization," IEEE Trans, commun., vol. 45, 

pp. 1497-1500, Dec. 1997. 
[17] Y. Li, N. Seshadri, and S. Ariyavisitakul, "Channel estimation for ofdm 

systems with transmitter diversity in mobile wireless channels," IEEE J. 

Select. Areas Commun, vol. 17, no. 3, pp. 461^47 1 , Mar. 1999. 
[18] E. Eweda, "Comparation of rls, lms and sign algorithms for tracking 

randomly time-varying channels," IEEE Trans. Signal Process., vol. 42, 

no. 11, pp. 2937-2944, Nov. 1994. 
[19] L. Lindbom, M. Sternad, and A. Ahlen, "Tracking of time-varying 

mobile radio channels- part i: the wiener lms algorithm," IEEE Trans. 

commun., vol. 49, no. 12, pp. 2207-2217, Dec. 2001. 
[20] M. K. Tsatsanis, G. B. Giannakis, and G. Zhou, "Estimation and 

equalization of fading channels with random coefficients," in Proc. of 

ICASSP'96, vol. 2, Atlanta, GA , USA, May. 1996, pp. 1093-1096. 
[21] G. C. Goodwin and K. S. Sin, Adaptive Filtering, Prediction, and 

Control. Englewood Cliffs, NJ: Prentice-Hall, 1984. 
[22] S. A. Hadei and P. Azmi, "A novel adaptive channel equalization method 

using variable step-size partial rank algorithm," in Proc. of AICT'10, 

Barcelona, Spain, May. 2010, pp. 201-206. 
[23] N. R. Yousef and A. H. Sayed, "Ability of adaptive filters to track carrier 

offsets and channel nonstationarities," IEEE Trans. Signal Process., 

vol. 50, no. 7, pp. 1533-1544, Jul. 2002. 



